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Abstract 

The  effects  of  irregular,  open  strnctnras  of  two  dimensional  aggregates  on  tha 
rates  of  diffusion  controlled  reactions  are  examined  by  computer  simulation. 
The  simulation  model  describes  the  reaction  between  two  species,  A  and  B,  where 
A  is  randomly  distributed  on  a  two  dimensional  square  lattice  and  diffuses  like 
a  random  walker  while  B  forms  a  fixed  aggregate  of  different  shapes  and  sixes. 
Both  unimolecular  (i.e.,  B  is  non-destructive  but  A  disappears  upon  reaction) 
and  bimolecular  (i.e.,  both  A  and  B  disappear  after  the  reaction)  mechanisms  are 
examined.  Different  sizes  of  DLA  aggregates  are  chosen  to  model  the  irregular 
aggregates  and  to  compare  the  reaction  rstes  with  the  circular,  compact 
aggregates  of  different  radii.  For  unimolecular  reaction,  with  the  same  total 
mass,  the  DLA  aggregate  shows  s  faster  reaction  rate  than  the  circular  one  for 
its  open  structure.  In  the  bimolecular  mechanism,  the  sizes  (or  radii)  of  the 
*ltreSate*  decrease  during  reaction  processes,  and  the  reaction  rate  is  slower 
compared  to  the  unimolecular  reaction.  This  shrinkage  effect  is  more  distinct 


for  DLA  i||r«|it«i  than  the  compact  aggregates.  It  la  foond  that  the  reaction 
ratea  of  DLA  aggregatea  are  compatible  with  that  of  circnlar  aggregatea  with  aa 
effective  radina.  The  relation  between  the  effective  radina  and  the  aaaaea  of 
the  DLA  aggregatea  can  be  described  through  a  power  law,  -  N0®.  where 

is  the  radios  of  circnlsr  aggregate  with  approximately  the  sane  reaction  rate  aa 
the  DLA  aggregate  of  total  aass,  Ne.  The  valnea  for  6  are  0.68  +  0.07  and 
0.79  +  0.08  for  nniaolecnlar  and  biaolecnlar  reactions,  respectively.  Froa  the 
effective  radios,  the  oniaolecolar  reaction  rates  can  be  calcolated  by  osing 
known  theories  of  diffosion  controlled  reaction. 


2 


I.  Introduction 


The  kinetic  stadia*  of  raaetioa  rata  prooasaaa  occuring  oa  solid 
surfaces1'1  have  bean  an  interesting  topic  for  different  fields  such  as 
reactions  on  catalyst,1-*  adsorption,  desorption  and  growth  processes®*1  on 
aetal  or  semiconductor  snrfaces.  Unlike  gas  phase  reactions,  where  the  reaction 
rates  sre  determined  by  the  collision  frequency,-  many  heterogeneons  reactions 
occnring  at  the  surface  of  solids  manifest  that  diffusion  is  an  important  step 
in  the  reactions  process. 1-1 

O 

Recent  studies  have  shown  that  many  dendritic  aggregates  formed  on  the 

surface  por sess  fractal  structure  and  obey  a  scaling  law,  N(R)  -  R®,  where  N(R) 

is  the  number  of  particles  within  the  radius  R,  and  D  is  the  fractal 

dimensionality.14  The  shape  and  the  fitted  fractal  dimensionality.  D,  of  these 

dendritic  aggregates  are  similar  to  DLA  (diffusion  limited  aggregate)  generated 

a 

by  computer  simulation  according  to  the  Vitten-Sander  model.  Because  many 
molecular  aggregates  involved  in  reaction  processes  possess  complex,  irregular 
structures,  it  is  difficult  to  model  these  systems  by  a  theoretical  approach. 
Therefore,  many  computer  simulation  studies4'®  have  been  carried  out  to  mimic 
possible  reaction  processes  occnring  on  the  two  dimensional  surface  with 
irregular  aggregates.  Silverberg  et.  a 1 . *  have  studied  the  effects  of 
adsorbate  interaction  and  aggregation  on  the  rates  of  chemical  reactions  on 
surfaces.  The  geometric  properties  (size  and  shape)  of  the  aggregates,  and  the 
lateral  interactions  are  considered  in  their  simulation.  In  their  model,  the 
migration  rate  of  reactant  is  dependent  on  the  assoc tsted  energy  changes  for 
each  step.  Therefore,  it  is  not  quite  clear  whether  the  process  is  'diffusion* 
controlled,  or  'reaction*  controlled. 
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In  the  present  work,  we  study  the  effect  of  diffusion  controlled  process  on 
the  surface  reactions  involving  aggregates  of  irregular  shapes.  In  this  case, 
the  contact  area  between  reactants  is  not  a  unique  factor  detenaining  the 
reaction  rate.  The  geometric  structure  would  play  an  important  role  for  the 
diffusion  process  snd  affect  the  reaction  rate.  In  this  study,  we  have  chosen 
the  DLA  to  model  the  dendritic  aggregates  of  irregular  shapes,  and  use  computer 
simulation  to  investigate  the  reactions  on  the  surface.  Te  assume  that  each 
particle  (atom  or  molecule)  is  a  hard  sphere,  and  no  long  range  interaction  is 
present.  When  two  similar  diffusive  particles  (A-A)  touch  each  other,  they 
repel,  and  reaction  happens  when  a  diffusive  psrticle  A  reaches  the  aggregate  B. 
Two  reaction  mechanisms  are  considered  in  the  simulation.  The  first  one  is 
unimolecular  process,  (e.g.,  catalytic  unimolecular  decomposition,  or 
isomerization),  in  which  the  aggregates  keep  the  same  shape  and  structure  during 
the  reaction.  The  second  one  is  bimolecular  process,  where  both  the  aggregate 
and  diffusive  particles  disappear  after  the  reaction.  Because  of  the 
non-uniform  shape  of  the  aggregates,  and  the  geometric  structure  of  the 
aggregates  play  an  important  role  in  determing  the  reaction  rate. 

The  concept  of  effective  radius,  refg,  is  used  to  compare  the  simulation 
results  of  the  DLA  clusters  to  circular,  compact  aggregates  where  many 
theoretical  studies  have  been  carried  out.*®  According  to  the  diffusion 
controlled  reaction  theories, *®  the  reaction  rate  would  increase  as  the 
aggregate  radius  increases.  Because  of  the  open  structure  of  the  DLA  aggregates 
(Fig.  1),  the  maximum  length,  rBIX#  from  a  particle  on  the  aggregate  to  the 
center  is  greater  than  the  radius  of  the  circulsr  aggregate  with  the  same  total 


mass.  However,  the  open  structure  of  DLA  aggregste  also  indicates  a  penetration 
effect  which  does  not  hsppen  for  the  circulsr  aggregates.  This  effect  would  let 
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the  probing  perticle  to  diffuse  among  the  open  spaces  between  the  dendritic 
branches  without  reacting.  Therefore,  the  effective  radius,  Ttff  would  be 
smaller  than  rBaz  hut  larger  than  the  corresponding  radius  of  a  circular 
aggregates  with  the  sane  total  mass. 

In  simulation  studies  of  the  growth  process  of  DLA  aggregate,11  an  active 

zone  was  found  where  most  new  particle  join  the  existing  cluster.  The  mean 

deposition  radius,  <r^>  shows  a  power  law  dependence  on  the  total  mass  of  the 

aggregate,  <r^>  ~  N*'’  .  The  value  of  J  for  DLA  aggregate  is  0.S84  +  0.02.  Our 

simulation  of  the  chemical  reaction  is  the  reverse  of  growth  process.  From  the 

simulation  results,  the  effective  radius  also  exhibits  a  scaling  law  behavior 

with  respect  to  the  total  mass  of  the  aggregate,  r  n  The  values  of  6 

ox  t  o 

for  unimolecular  and  bimolecular  mechanisms  are  0.68  +  0.07  and  0.79  +  0.08, 
respectively. 

II.  Theoretical  Background  of  the  Diffusion  Controlled  Reactions: 

Consider  reactions  between  reactants  A  and  B  with  the  following  reaction 
equations : 

A  +  B  - >  B  +  C  (la) 

A  +  B - >  C  (lb) 

where  both  species  (A  and  B)  are  randomly  distributed  and  the  reaction  product  C 
has  no  influence  on  the  further  reaction  processes.  Both  of  two  cases,  (la)  and 
(lb),  are  diffusion  controlled  processes.1®  The  first  case  is  like  a  catalytic 
reaction,  where  A  particles  are  mobile  and  destructible  while  B  is  the  fixed 


aggregate  (e.g.,  a  catalyst)  and  is  indestrcntible.  An  example  of  this  conld  be 
a  catalytic  nniaolecnlar  chemical  process  (e.g..  nniaolecnlar  decoaposition  or 
isomerization).  Assuming  A  particles  are  undergoing  an  isotropic  diffusion  in 
tvo  dimensional  lattice  like  a  random  walker,  after  several  steps,  an  A  particle 
reaches  the  aggregate  B.  and  reacts  (i.e.,  disappears)  in  the  simulation.  By 
counting  all  the  survival  A  particles  at  different  time  t.  the  reaction  rate  can 
be  calculated  by  the  following  equation 

_ JCA<«> _ 

dt 


Rate 


(2) 


Cg(t)  =  constant 

7.  5 

The  second  case  is  a  bimolecular  surface  reaction:  '■*  e.g..  0  +  CO - >  CO^ .  At 

the  beginning  of  the  reaction,  one  of  the  species,  B,  forms  an  aggregate  on  the 
surface.  The  other  reactant  A  covers  the  surface  randomly.  These  A  particles, 
then,  start  to  diffuse  on  the  surface  like  random  walkers.  Once  they  reach  the 
aggregate  B,  both  A  and  B  react  and  disappear  in  the  simulation.  Therefore,  the 
size  of  aggregate  B  will  shrink  as  the  reaction  goes  on.  By  recording  the 
changes  of  the  concentration  of  A  (or  B)  with  time,  the  reaction  rate  law  can  be 
calculated  as 


Rate 


dCA(t)  dCB(t) 

dt  dt 


(3) 


To  simplify  the  problem,  we  take  the  the  center  of  the  aggregate  B  as  the 
origin  of  the  coordinates  of  the  A  particles.  Assuming  B  particles  are  immobile 
and  indestructible  sinks  and  A  particles  undergo  isotropic  two  dimensional 
diffusion  with  a  diffusion  constant  D,  then  the  spatial  distribution  of  the  A 
particles  around  a  sink  (aggregate  B)  can  be  described  in  terms  of  a  space  and 
time  dependent  function  pA(r,t),  normalized  with  respect  to  the  mean 
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concentration  C^ft).  pfe(r,t)  satisfies  the  diffusion  equation: 
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3  PA(r,t) 

d  t 


D  2pA(r,t) 


(4) 


with  appropriate  initial  and  boundary  conditions. 

If  the  system  is  spherically  symmetric,  the  reactants  A  and  B  are  assumed 
spherical  with  radii  r^  and  rg,  respectively.  A  reaction  takes  place  if  two 
particles  approach  one  another  within  the  contact  distance  r^g,  «  r  +  tg. 

By  using  Laplace  transformation,  the  partial  differential  equation  (Eq.  4)  can 
be  solved  from  the  following  boundary  conditions: 

p(oO,t)=l.  (5a) 

p(rAB't)=0*  for  t  >  0  (5b) 

p(r,0)*l.  for  r  >  r^  (5c) 

The  reaction  rate  can  be  obtained  from  the  diffusion  flow  of  A  particles 
into  the  boundary  of  the  sinks  B  with  concentration  Cg, 


4CA(t) 


dt 


—  -  -  k(t)  Cg  CA(t) 


(6) 


where  the  rate  constant  k(t)  can  be  described  as*®^ 


d  p(r,t) 

k  ( t )  =  2nDrAg  ( - ) 

3  x 


2-dim 


r-r 


AB 


(7) 


which  decreases  as  time  increases  due  to  the  concentration  depletion  on  the 
boundary.  The  reaction  rate  k(t)  is  also  dependent  on  the  contact  distance  r^g. 
The  larger  the  radius  (or  the  contact  surface  area)  is,  the  faster  the  reaction 
rate  would  be. 
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1  3 

Owen1*  bis  derived  ia  analytic  solution  for  the  2-dimensional  esse  in 
cylindrical  coordinates  and  gave  a  closed  fora  of  k(t)  as  a  function  of  r^,  and 

D.  In  a  good  approxiaation,  a  simple  interpolation  foranla  was  found  for 
k(t)>10(d) 

k(t)  -  nD  (a(t )  +  2^  („Dt)-1/2)  (8) 

o(t)  =  Min  (0.6,4/ln  (4Dt  /  (nr^2) ) ) 

However,  for  the  aggregates  of  irregular  shape,  there  is  no  analytic  solution. 
Thus,  the  concept  of  effective  rsdius,  r  which  plays  the  role  of  r^  in  the 
above  equation  is  proposed  to  describe  the  reaction  rates  of  DLA  aggregates. 
Once  the  effective  radius  is  obtained,  the  unimolecular  reaction  rate  can  be 
calculated  through  the  existing  solutions  for  two  dimensional  disks ,10^d ) 

III.  Simulation  of  Unimolecular  and  Bimolecular  Reaction 

The  simulation  is  carried  out  on  the  two  dimensional  square  lattice  within 
a  circular  boundary  of  300  lattice  unit  in  diameter.  The  procedures  of 
simulation  are  as  follows:  First,  a  aggregate  B,  which  can  be  any  shapes  and 
sizes,  is  generated  in  the  center  of  the  two  dimensional  lattice.  At  time  equal 
to  zero,  the  A  particles  with  a  given  concentration  are  randomly  distributed  on 
the  same  lattice  where  B  already  exists.  Then,  let  A  particles  diffuse  randomly 
like  a  random  walker.  Once  particles  A  reach  the  aggregate  B,  they  react. 
Either  A  alone  or  both  A  and  B  disappear  depending  on  the  mechanisms  we  are 
considering.  During  the  reaction  period,  we  keep  tracing  the  positions  of  both 


species  end  record  the  ntusber  of  sll  the  survival  (unreacted)  particles  (i.e., 
concentration)  to  calculate  the  reaction  rate. 

Doe  to  the  tine  limitation  of  coapnter  simulations,  the  largest  radius  of 
the  aggregates  used  for  this  simulation  is  60,  which  is  reasonable  compared  to 
the  radius  of  the  whole  system,  ISO.  The  total  number  of  steps  for  the 
simulation  is  2000,  and  the  average  distance  travels  by  a  randomly  diffusing 
particle  is  2000*^*  **  4S  lattice  units.  Therefore,  the  results  of  simulation 
would  not  be  perturbed  by  the  finite  sire  effect.  Each  initial  condition  (e.g., 
the  aggregate  shape,  and  size)  is  run  three  times  and  the  results  are  averaged 
to  reduce  the  fluctuations  in  the  reaction  rate  calculation.  When  a  particle 
reaches  the  circular  boundary,  it  would  be  deleted  and  a  new  particle  is 
generated  on  the  opposite  site  of  the  boundary.  The  decay  of  the  reactant 
concentration,  CA(t).  is  normalized  to  its  initial  concentration,  C^(t  ■  0)  to 
compare  the  reaction  rate  for  different  initial  conditions. 

The  simulations  sre  csrried  out  with  two  kinds  of  aggregates  (Fig.  2). 
First,  circular,  compact  aggregates  with  different  radii,  which  are  equivalent 
to  the  contact  distance  r^g  in  equation  7,  are  used.  The  simulation  results  are 
then  used  as  the  reference  curves  to  determine  'reff'  for  the  aggregates  of 
irregular  shapes.  The  aggregates  generated  by  the  diffusion  limited  aggregation 
process  are  chosen  to  model  the  irregular  aggregates  in  real  systems.  The  sizes 
of  DLA  aggregate  are  varied  for  comparison  to  the  reaction  rates  of  the 
circular,  compact  aggregate.  The  effective  radius,  reff,  for  esch  DLA  aggregate 
(with  different  total  mass)  is  determined  by  comparing  its  reaction  rates  to  a 
series  of  compact  circular  aggregates  to  find  which  one  will  give  the  same 
reaction  rate.  By  plotting  log  NQ(reff)  vs.  log  (reff),  a  straight  line  is 
obtained  and  indicates  a  scaling  law  behavior.  The  simulations  are  carried  out 


for  both  the  unimolecular  and  bimolecular  processes,  and  the  sealing  exponents 
are  determined  respectively. 

IV.  Results 

The  simulation  consists  of  calculating  C^(t)  at  different  tiaie  t,  for 
different  shapes  of  aggregates  (i.e..  circular  or  DLA)  and  for  different  sisea 
of  aggregates.  The  simulation  is  performed  on  IBM  3090  computer,  and  follows 
the  procedure  described  in  section  III. 

A:  Effect  of  The  Aggregate  Shape: 

The  decay  of  the  concentration  of  A  at  different  time  with  differen  sites 
(or  radii)  of  the  aggregates  are  shown  in  part  A  of  Fig.  3.  and  4  for  circnlar 
aggregates,  and  part  B  of  Fig.  3  and  4  for  DLA  aggregates.  In  part  A,  the  total 
mass  of  DLA  aggregates  are  2000,  1800,  1300.  1300,  1000,  800,  and  300.  In  part 
B.  the  radii  (total  mass)  of  the  circle  are  60  (11310),  33  (9303).  30  (7834),  43 
(6362),  40  (3027),  33  (3848),  30  (2827),  and  20  (1237)  lattice  units.  The 

concentration  of  the  diffusive  psrticle  A  is  0.06.  For  the  same  total  mass,  the 
DLA  aggregate  reacts  faster  then  the  circular  one.  In  other  words,  it  has 
larger  effective  radius.  r#ff.  The  plot  of  r#ff  vs.  is  given  in  Fig.  3. 

The  slope  of  straight  lines,  6,  is  equal  to  0.68  +  0.07  and  0.79  ♦  0.08  for 
unimolecular  and  bimolecular  reactions,  respectively. 

B:  Effect  of  the  Reaction  Mechanism,  Unimolecular  vs.  Bimolecular: 


Ia  Fig.  I  lid  4.  tk«  reaction  Mckuiiai  arc  MiaalMiUr  a *4  klatUnUr, 
respectively.  For  tk«  biswlecular  aackaaiia,  tka  aiaa  (or  radtaa)  of  tk« 
aggragata  decreases  as  time  iaaraasas.  Tkti  would  alow  4owa  tka  raactiow  ta 
aoafariic.  to  tka  uaimolecalar  mechanism.  Tka  skriakaga  affaet  ia  aor a  giatiaat 
for  DLA  aggregates  tkaa  cirealar  oaaa  4aa  to  its  op# a  atraetara. 

V.  Discassioa 

Froa  abova  plots,  wa  caa  saa  tkat  tka  akapa  of  tiaai  4apaa4oat  eoaceat  rat  ioa 
faactioa  of  raactaat  A  for  DLA  aggragataa  la  aiailar  to  tkat  of  coapact. 
cirealar  aggragataa.  Bowever,  tka  DLA  aggragataa  skow  fastsr  raactioa  rata  tkaa 
tka  circular  oaa  witk  tka  aaaa  aggragata  sita  (total  boss).  Tkis  difference  la 
mostly  4aa  to  tka  geometric  atraetara  of  tka  aggragataa  wklak  will  ckaaga  tka 
boundary  coaditioas  of  tka  diffarcatial  equation.  Ig.  4  a ad  3.  Basically,  tka 
raactioa  kinetic  for  the  dendritic  DLA  aggragata  caa  ka  described  by  tka  aaaia 
diffareatial  aquation.  Eg.  4.  bat  with  different  boundary  conditions.  As  shown 
ia  Fig.  2B,  for  tbs  compact,  circular  aggregate,  as  tka  diffusing  raactaat  A 
raackas  tka  circular  boundary  of  aggragata  B,  the  raactioa  occurs  immediately. 
Tkis  implies  tka  same  boundary  coaditioas  aa  shown  in  Eg.  3.  Bowavar.  for  tka 
irregular  shape.  DLA  aggragata.  (Fig.  2A),  when  raactaat  A  raackas  tka  region  1. 
it  caa  altkar  react  whaa  reach  aggregate  B  or  still  ka  alive  and  diffuses  among 
tka  dendritic  braaebas.  Therefore,  the  DLA  aggragata  possesses  complicated 
boundary  conditions  for  the  differential  equation  and  does  not  have  analytic 
solut ioa. 

For  tka  DLA  aggregate,  because  of  its  open  struetura,  the  mar  imam  'length' 
of  brauckas  which  determines  the  radius  of  the  boundary  between  regioa  I  and  II 


(Fig .  2A>  is  ii 


lirculi 


with  tk( 


im«  Mout  of  total  uii.  Matheaatically,  this  behavior  can  bo  found  from  tba 
aoaliag  lav  for  fractal  structures,  N(l)  -  1®.^  Tko  acaliag  law  caa  bo 
rewritten  aa 


raax  '  <»> 

For  the  fractal  structures,  it  is  kaova  tkat  d  2.  D.  Therefore,  tke  aaxiaua 

leagtk  r>M<  is  larger  tkaa  tke  corresponding  radios  of  coapact  disk.  Siace  tko 

value  of  ra  is  related  to  the  radios  of  tke  boundary  betweea  region  I  and  II, 

which  is  equivalent  to  tke  contact  distance  r^g  as  ia  eqaatioa  8.  the  larger  tke 

aaxiaua  lengtk  r^^  is,  the  faster  the  reaction  rate  would  be.  However,  due  to 

tke  open  structure  of  tke  OLA  pattern,  wkea  tke  diffusive  particle  A  reaches  tke 

distance  which  is  equal  to  tke  aaxiaaa  leagtk  of  tke  branckes,  it  nay  not 


:ily  react  at  tkat  timt.  Tke  particle  caa  still  diffuse  but  witk  tke 


sarvivial  probability  decreasing  with  tine.  This  indicates  that  the  real 
aaaber  of  reacted  particles  A  is  less  tkaa  tkat  witk  tke  circular  boundary  for 
tke  saae  boundary  radius.  Therefore,  these  two  conflicting  factors:  larger 

aaxiaua  lcagt’  and  penetratiag  effect,  would  coapete  each  other,  and  the 
siaulatioa  results  show  tkat  the  iacreaaa  ia  tke  aaxiaua  length,  r^  ,  is  aore 
iaportaataat  than  tke  effect  of  tke  penetratioa  aaoag  branckes  for  tke  OLA 
aggregetes.  By  coapariag  tke  reaction  rate  of  DLA  aggregate  witk  given  total 


aass  vitk  reaction  rates  of  a  series  of  circular  aggregates  witk  differeat 


effective  radius,  r 


cff,  for  OLA  aggregate  (which  is  equivalent  to 


tke  contact  radius  r^g  is  Eq.  I)  can  be  obtained.  rcff  *•  deterained  as  tke 
length  which  is  eqael  to  the  radius  of  circular  aggregate  witk  tke  approxiaate 

the  saae  reectioa  rate  as  DLA  aggregate  of  gives  total  aass,  N  Because  of  tke 

1/2 

dendritic,  open  streture  of  DLA  aggregate,  is  larger  tkaa  N#  (i.e.,  tke 
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radius  of  disk).  However,  reff  is  smeller  then  raaz  because  tbe  diffusing 
reactant  A  nay  diffuse  anong  the  dendritic  branches,  and  penetrates  closer  to 
the  center  of  aggregate. 

The  effect  of  reaction  mechanisms,  uniaolecular  vs.  bimolecular,  is  also 
examined  in  the  simulation.  For  the  bimonlecular  case,  the  aggregate  sites  will 
disappear  when  reached  by  the  diffusive  A  particles,  and  the  effective  radius 
r9f(  of  the  aggregate  will  shrink  as  reaction  goes  on.  Since  the  time  dependent 
reaction  rate  constant,  k(t),  depends  on  r^  which  is  decreasing  with  time  for 
bimolecular  mechanism  and  is  a  constant  for  nnimolecular  case,  the  bimolecular 
reaction  rate  is  slower  than  the  nnimolecular  mechanism.  Because  of  the  open 
structure  of  dendritic,  DLA  aggregate,  for  the  same  amount  of  particle  B 
(aggregate)  consumed,  the  change  in  the  maximum  length,  r  |s  larger  for  the 
DLA  than  the  circular  aggregate.  Therefore,  the  shrinkage  effect  is  a  more 
dominant  for  the  DLA  than  the  circular  aggregates,  and  for  the  aggregate  with 
smaller  sizes  than  the  one  with  larger  sizes.  Due  to  the  complexity  of  the 
shrinkage  effect  on  the  geometric  boundary  of  the  aggregates  for  bimolecular 
mechanism,  and  the  random  positions  of  the  aggregate  sites  being  reacted,  there 
is  no  analytic  solution  available  even  for  the  compact,  circular  aggregates. 
The  results  of  reff  for  bimolecular  mechanism  can  only  describe  qualitative 
behavior,  and  can  not  give  actual  reaction  rate. 

o 

In  the  simulation  of  growth  kinetic  of  DLA  aggregate,  the  growth  can  be 
described  in  terms  of  an  active  zone  defined  as  the  region  where  new  particle 
join  the  existing  cluster.  This  zone  is  characterized  that  the  Nth  particle  is 
deposited  to  the  aggregate  with  (N  -  1)  particles  at  a  distance  r^  froB  the 
center  of  mass.  The  mean  deposition  radius  <r^>  shows  a  power  law  scaling 
behavior  to  the  total  mass  of  the  aggregate. 
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The  vslue  of  O  for  DLA  aggregate  ia  0.584  +  0.02.®  Is  this  paper,  we  study  the 
diffusion  controlled  reaction  which  is  the  reverse  of  the  growth  process.  If  we 
plot  the  versus  NQ  in  the  log  log  scale,  a  straight  line  can  be  obtained  es 
shown  in  Fig.  5  for  both  unimolecular  and  bimolecular  mechanism.  The  slope,  6, 
is  equal  to  0.d8  +  0.07  and  0.79  +  0.08  for  uni-  and  bi-nolecular  aechanisa, 
respectively.  The  difference  in  the  values  of  J  and  6  (for  uniaolecular 
aechanisa)  between  the  growth  and  the  reaction  process  could  be  due  to  the 
initial  distribution  of  the  diffusive  particles  A.  For  the  growth  process,  all 
the  probing  particles  are  staying  at  far  away  from  the  aggregate.  However,  in 
our  simulation,  the  probing  particles  A  are  randomly  distributed  on  the  whole 
lattice  initially.  This  would  sllow  some  particles  staying  among  the  branches 
of  the  aggregates  in  the  beginning  of  the  reaction.  These  particles  conld  react 
with  the  sites  close  to  the  center  of  aggregate  without  feeling  the  acreening 
effect  from  the  outer  branches.  The  meaning  of  6  for  bimolecular  mechanism  is 
not  very  clear,  and  it  is  difficult  to  compare  the  values  of  and  6. 

For  the  bimolecular  mechanism,  the  reaction  consumes  the  same  amount  of 
reactants  A  and  B.  If  the  number  of  reactant  A  reacted  is  compatible  with  the 
total  number  of  aggregate  B,  the  relative  change  of  r^^  would  be  larger  during 
the  reaction  period,  and  the  shrinkage  effect  is  more  dominant  in  the  smaller 
aggregate.  These  factors  would  lead  to  a  smaller  value  of  8  compared  to  the 
unimolecular  case.  However,  the  total  amount  of  reactant  A  (initial 
concentration)  would  also  affect  the  exponent  &.  Because  of  the  complexicity  of 
the  shrinksge  effect,  there  is  no  theoretical  equation  to  predict  the  reaction 
rate  even  for  the  circular  boundary.  The  comparison  between  two  kinds  of 


aggregates  gives  only  the  qualitative  results,  snd  there  is  no  direct 


correlation  between  the  values  of  the  6  exponent  for  uni-  and  bi-aolecnlar 
■echaniaa. 

-  < 

In  the  earlier  aianlation  of  Silverberg  et.  al.,  they  have  etodied  the 

effect  of  the  geometric  propertiee  of  aggregates  to  the  reaction  kinetica.  The 
reaction  rate  is  dependent  on  the  nnaber  of  surface  sites  which  are  available 
for  the  diffusive  particles  to  react.  However,  for  the  dendritic  DLA 
aggregates,  beesnse  of  the  strong  screening  effect,  even  though  the  surface  area 
is  linear  to  the  total  aass  of  the  aggregate, only  those  sites  which  are 
located  at  outer  branches  of  the  aggregate  are  useful  for  the  reaction  to  occur. 
Therefore,  the  diffusing  reactant  A  slight  feel  different  screen  effect  of  the 
irregular  surface  (boundary)  for  the  'diffusion'  controlled  process  and  the 
'reaction'  controll  process.  In  the  present  study,  the  'diffusion'  controlled 
process  is  investigated,  and  we  have  concentrated  on  how  the  radii  of  these 
'active  zone’  are  varied  during  the  reaction  with  the  size  of  aggregates.  Vith 
the  help  of  known  results  of  diffusion  controlled  reaction  on  compact,  circular 
disk  (Eq.  8).  the  inforastion  of  reaction  rate  can  be  obtained  f row  a  single 
parameter,  *eff,  for  the  irregular  shape  aggregatea. 
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Figure  Captions: 


Figure  1: 

The  aggregate  gernerated  by  the  diffusion  United  aggregation  process  as 
described  in  reference  9.  The  total  mass  of  the  aggregate  is  12.000  lattice 
points. 

Figure  2: 

The  slape  of  aggregate  used  for  the  simulation.  (A)  DLA  aggregate  with 
2000  lattice  point.  (B)  Circular  aggregate  with  radii  equal  to  50. 

Figure  3: 

Plot  of  the  reactant  A  concentration  C^(t)  at  different  time  t.  C A(t)  is 
normalized  to  the  initial  concentration  C^(t>0) .  The  simulation  is  carried  out 
for  unimolecular  mechanism  for  the  DLA  aggregate  (part  A)  and  circular  aggregate 

(part  B) .  In  part  A,  the  total  mass  of  DLA  aggregates  are. -  2000, 

- 1800. -  1500. -  1300, -  1000. - 800.  and 

-  500,  respectively.  In  part  B,  the  radii  (total  mass)  of  circular 

aggregates  are  -  60  (11310). -  55  (9503), -  50  (7854), 

- 45  (  6362). -  40  (  5027), -  35  (  3848), -  30  (  2827).  and 

■  20  (1257),  respectively.  For  the  same  total  mass,  the  DLA  aggregates 
react  much  faster  than  the  circular  one  i.e.,  DLA  aggregate  has  larger  effective 
radius . 


Figure  4: 


Plot  of  the  reactant  A  concentration  C^(t)  at  different  tiae.  C^(t )  ia 
normalized  to  the  initial  concentration  C^(t«0).  The  aimnlation  ia  carried  ont 
for  bimolecnlar  mechanism  for  the  DLA  aggregate  (part  A)  and  circnlar  aggregate 
(part  B) .  In  part  A,  the  total  mas a  of  DLA  aggregatea  are  -  2000. 

-  1800. -  1500. - 1300, -  1000. -  800,  amd 

500,  respectively.  In  part  B,  the  radii  (total  mass)  of  circnlar 

aggregates  are  -  60  (11310). -  55  (9503). - 50  (7854). 

- 45  (6362), - 40  (  5027). - 35  (  3848),  -  30  (  2827).  and 

* - —  20  (1257).  respectively.  For  the  binolecular  mechanism,  the  size  (or 

radios)  of  the  aggregate  decreases  as  time  increases  which  would  alow  down  the 
reaction  in  comparison  to  the  unimolecular  mechanism.  The  shrinkage  effect  ia 
more  dominant  for  DLA  aggregate  than  the  circnlar  one  dne  to  its  open  structure. 

Figure  5: 

The  log-log  plot  of  with  N0#  where  *tff  ia  the  radius  of  circular 

aggregate  with  the  same  reaction  rate  as  the  DLA  aggregate  of  total  mass. 
Plots  A  and  B  are  the  results  of  unimolecular  and  bimolecnlar  mechanism 
respectively. 
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